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We construct a consistent action for a massive spinning quark on the end of a QCD string that 
leads to pure Thomas precession of the quark's spin. The string action is modified by the addition 
of Grassmann degrees of freedom to the string such that the equations of motion for the quark spin 
follow from boundary conditions, just as do those for the quark's position. 



I. INTRODUCTION 

A consistent description of spin within a QCD string 
theory has been sought for many years. The addition of 
dynamical spin to the bosonic string led to the develop- 
ment of supersymmetry and superstring theory 0] . Such 
theories are more realistic as unified theories of elemen- 
tary particle physics than as phenomenological descrip- 
tions of hadronic states. 

A more realistic description of hadronic states involves 
the replacement of the free end of the dual resonance 
string by the addition of a massive point quark to the 
end of the string. In 1977 Ida analyzed the motion of 
a spinless massive quark on the end of a bosonic string. 
The relativistic flux tube model , derived from different 
assumptions, is mathematically equivalent to a bosonic 
string with a spinless quark end and produces realistic 
meson spectra on average, but there is no place for quark 
spin in this model. In this paper we make a modification 
of the bosonic string plus bosonic quark model to intro- 
duce quark spin. 

Our clue to constructing a consistent action comes 
from the suggestion of Buchmuller [j] that the spin of 
the quark should undergo pure Thomas precession be- 
cause the quark sees a purely chromoelectric field in its 
rest frame. This seems to be supported by experimental 
data 0, |(J and is in agreement with QCD 0, ■ 

We begin in Sec.[n]by discussing the treatment of spin 
in pseudoclassical language. We show how to construct 
actions for a free fermion as well as a fermion with back- 
ground scalar and vector potentials. We analyze the case 
of a scalar potential in detail and show how the Thomas 
precession manifests itself in this language. 

In Sec. IIIII we show in detail that the Fermi- Walker 
transport of the spin vector, which is the equation of 
motion of the spin vector for a particle in a scalar po- 
tential, leads to Thomas precession of the spin in its rest 



frame. 

In Sec. II VI we use the example of a spinless quark 
coupled to the end of a string to argue for the form of 
the action for a spinning particle coupled to a modified 
Polyakov string action. The key idea is to obtain the 
equations of motion of the spin of the quark from bound- 
ary conditions, just as the equations of motion of the 
quark's position arise from boundary conditions. To this 
end, we introduce new Grassmann-valued fields on the 
string worldshoet. 

In Sec. El we use the consistency of the equations of 
motion of the quark and the requirement of Thomas pre- 
cession to fix the parameters in the string action. The 
result is that the only modification of a free spinning 
quark plus free bosonic string action is the replacement 
of the bosonic string position variable by the string posi- 
tion variable plus a term bilinear in worldshcct fcrmionic 
variables. 

In Sec. lVII we explore the fermionic gauge invariance of 
our string action. In the phenomenologically interesting 
case, we find that the worldsheet fermionic variables are 
pure gauge degrees of freedom. 

We find the momentum and angular momentum from 
Noether's theorem in Sec. IVIII These conserved quanti- 
ties are the usual starting point for the numerical quan- 
tization of the relativistic flux tube model. Finally, we 
conclude in Sec. IVIIII 



II. SPIN IN PSEUDOCLASSICAL MECHANICS 

We choose to work within the framework of pseudo- 
classical mechanics because the formalism is elegant 
as well as physically transparent; the transition from 
pseudoclassical to quantum mechanics is immediate. In 
this section we construct actions that produce the Dirac 
equation, both free and in background potentials, as an 
equation of motion and we show how the Thomas preces- 
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sion in a scalar potential manifests itself in this language. 
The main disadvantage is that it requires some familiarity 
with the technical details of Dirac's constrained Hamil- 
tonian mechanics as well as classical mechanics 

with Grassmann variables 0, |nj . 

The easiest way to construct pseudoclassical actions 
for fermions is to consider the Dirac equation as a phase- 
space constraint and to construct consistent actions that 
yield this constraint. The first actions of this type were 
found by Berezin and Marinov [l^. Barducci, Casal- 
buoni, and Lusanna Il3l . and Brink, Deser, Zumino, Di 
Vecchia, and Howe 14] . To represent the spin degrees 
of freedom of a fermion, a set of five Grassmann coordi- 
nates, £ M and £5, are introduced. Upon quantization, the 
Grassmann coordinates will become generators of a Clif- 
ford algebra and can be identified with Dirac's gamma 
matrices. The kinetic piece of the action for the Grass- 
mann variables 



/ dr l - + f 6 6 



■Vku,.,,, = / ./"- U„t'' + • (2.1) 

leads to the canonical second-class constraints 



0, 



X5 = 7T S 



2 & 



0. 



(2.2) 



Here we use Dirac's wavy equal sign notation [Tol ITU 
for "weak equality," which reminds us that the equalities 
cannot be taken before Poisson brackets are calculated. 
We denote the canonical momenta to £^ and £5, defined 
to be the derivative of the Lagrangian from the right 
with respect to the velocities £ M and £5 respectively, by 
and 775 . With this convention, we obtain the following 
Poisson brackets 



{£5,^5} = {7r 5 ,£ 5 } = 1, 



(2.3) 
(2.4) 



with all others being zero. Our conventions for pseudo- 
classical mechanics are given in Appendix 1X1 

The weak equalities in Eq. I|2.2[) can be replaced by 
strong ones if we introduce the Dirac brackets [To| • From 
the definition in Appendix [B] and the Poisson brackets 
above, we find 



&}£> = -iVlfl 

{^,^}d = 0, 

{£5,&}£> = ~i, 



(2.5) 
(2.6) 
(2.7) 



where rj^ is the metric. Our convention is r/^ — 
diag(-l,+l,+l,+l). 

The meaning of the Grassmann numbers becomes clear 
upon quantization. When we make the replacement of ih 
times Dirac brackets by anticommutators, we find that 
the quantum operators £ M and £5 obey a Clifford algebra, 



£/j£^ ~\~ £^£^1 
&6 



0, 

h 
2' 



(2.8) 



From these anticommutation relations, we see that the 
operators £^ and £5 can be represented as gamma matri- 
ces, 



P I* 

£5 = Y2 75 - 

The free Dirac equation is proportional to 

= (p^ + m^) |V) =0. 

Thus, we should introduce the constraint 
4> = Put" + m£ 5 fts 



(2.9) 
(2.10) 

(2.11) 
(2.12) 



into our action. This constraint does not have vanishing 
Dirac bracket with itself, but yields the Klein-Gordon 
operator: 



K = ^{P^ +m^,p„e +m^ 5 } } 
= \{p 2 +m 2 )~0. 



(2.13) 



In order to be able to impose the constraint <fi as in 
Eq. (|2.11l) . cf> and any constraints, such as K, arising from 
it must be first-class, which means the Dirac brackets of 
any pair of them yields a combination of other first-class 
constraints. In order for the set of constraints to close 
under Dirac brackets, this last constraint must have van- 
ishing Dirac bracket with <j>. This is guaranteed by the 
(graded) Jacobi identity, 

{</>, K} D = {0, <j>} D } D = . (2.14) 

The dynamics of this system are given by the free ac- 
tion, plus these constraints put in with Lagrange multi- 
pliers A, and e: 



S = J dr[ Pll x* + + 

+ i-fat" + m&) - 4(p 2 + m 2 )j (2.15) 
m 2 J 

We may eliminate p from Eq. I|2.15l) . by using its 
(purely algebraic) equation of motion. Similarly, we may 
then eliminate e from the intermediate action to find the 
action given by Berezin and Marinov [12J 



S = 



(It 



+ iA(£ 5 + u ■ £) 



(2.16) 



where we have used the usual notation for the four- 
velocity, u — x/\J—x 2 . 
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The Dirac equation in a background scalar field ip and 
vector field is obtained from the free equation by min- 
imal substitution for and the addition of <p to the 
mass: 



<i> = (p- A) fi e + (m + <p)& «0 



(2.17) 



We wish to use this as a constraint to construct an 
action in the same manner. We must again consider that 
the constraint <f> have vanishing Dirac bracket with itself. 
We find 



K 



~{<M}d = ~(p- A) 2 + i(m + (^) 2 
-k^^+^^wO. (2.18) 



where i* 1 ^ = d^A^ — d u A^. The Jacobi identity again 
insures that there are no further constraints. 

As before, we implement these constraints by use of 
Lagrange multipliers, a commuting one, e, and an anti- 
commuting one, A, 

S = Jdr [p^ + L^ + frfa+i^-eK]. (2.19) 

We note that the action for a spinless particle can be 



obtained by taking the spin variables to zero: £ p 



0, 



6->o. 

Because we are interested only in the Thomas preces- 
sion here, from now on we consider the action with a 
scalar potential only, so we set = 0. Eliminating first 
Pft, and then e, in the action Eq. i|2.19[l with A^ = 0, we 
find 



S = 



(It 



ip- 



to + <p 



(^ + &&)+iA(e B +u-0 



(2.20) 



This action is the same as the one analyzed by Marte- 
myanov and Shchepkin We note that the Thomas- 
Bargmann-Michel-Telegdi equations of motion 061 IVA 
U8| for the spin can be found from an analysis |l2L Il5l 
El US E3 of the action $FWi with <p = and A^ ^ 0. 



The equations of motion from the action 
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F, 



Xu, 



where 



F„ 



ip- 



-9, 

dL 

dxf* 

(to + (p)u^ 
. A 



TO -f 
TO + ip ' 



'6s£ • dip 
m + ip 



are 
(2.21) 

(2.22) 
(2.23) 



(2.24) 



(2.25) 



and we have used the convenient notation 



(2.26) 



for the projection operator perpendicular to the four- 
velocity. 

In order to clarify the algebra in the rest of this section, 
we follow Martemyanov and Shchepkin |Tfij and work to 
lowest order in the fermionic variables. In this approxi- 
mation we have 



P F v 

TO + ip ' 



(2.27) 



u ■ F 

A I «,i (2.28) 



m + (p 

& = A ■ u-£. 

TO + ip 



(2.29) 



The momentum and angular momentum of the sys- 
tem can be found by Noether's theorem. We make an 
infinitesimal Poincare transformation of the variables 



Sx" = a»+uj^ v x v , 

and extract the conserved quantities from 
(dL 



(2.30) 



53 - HiHWf^ 

a»A Pfi + -u^AJ^, (2.31) 



where d /d^ denotes the derivative acting from the 
right, and A denotes the difference in values between 
final and initial times. In Eq. (|2.31(l we have also used 
the equations of motion. 

We find that the total angular momentum is a sum of 
orbital and spin pieces 



Jau Luy Sn U %\uPv\ ^SWSf ■ 



(2.32) 



The total angular momentum, as well as each piece sep- 
arately, obeys the Dirac brackets relation 

~\~ ~l~ Vfj'P' Jvv f • (2.33) 



The Pauli-Lubanski vector, 



(2.34) 



represents the spin of the particle and is purely spatial 
in the rest frame of the particle; 



u ■ s = 0. 



(2.35) 



We use the convention that £0123 = +1. Using the iden- 
tity 



(2.36) 
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we may revert Eq. I|2.34[l to find 

iff" = - e a ^ s u jS s + i{u a f - u P C)(u ■ 0- (2-37) 
Using Eq. (|2.37l) , we find the rate of change of s M 

+ ie^ af3 n u i a ^. (2.38) 

We observe that the equation of motion for f must have 
the form 

^ = -^{u ■ £) + ^(anything), (2.39) 

in order for the Pauli-Lubanski vector to be Fermi- Walker 
transported along the worldline of the particle. That is, 
for to obey 



L u v s v = (u^u v - u IJt u v )s v . 



(2.40) 



Equation (|2.40(1 is the condition that there is no torque 
on the spin. The spin thus undergoes Thomas precession, 
as we will see in the next section. 



The equation of motion satisfied by the rest-frame spin 
vector is 



= — (As) = As + As 
dt dt y ' 

= A A -1 s + AftA -1 s - 
The rotation matrix l|3.2|l is 



(3.6) 



W v = V K ■ ' ■ ) ■ (3-7) 
v vv — vv 



Simplifying the right hand side of Eq. 13.6(1 , we find 
ds 1-1 f 



dt » 2 \0 vv- vv 



s . 



(3.8) 



Since the rest frame spin, so, has no time component, we 
have 

ds Q 7-1 7 - 1 

— = — — (vv - vv) • s = 5— (v x v) x s . 

at v z v z 

(3-9) 

The acceleration of a particle in uniform circular motion 
with angular velocity ui is 



III. THOMAS PRECESSION 

In this section we demonstrate that a vector that un- 
dergoes Fermi- Walker transport in a circular orbit will 
precess in its rest frame at the Thomas frequency. 

The spin vector of a gyroscope moved along a space- 
time path (r) in the absence of net torque undergoes 
Fermi- Walker transport. We take laboratory time to be 
the worldline parameter; r = t. The rate of change of its 
spin vector then is 



ds^ 1 
~~dt 



with 



u^u v — wu Vl 



(3.1) 



(3.2) 



where the is the four velocity tangent to x^(t) and 
dot means derivative with respect to t. 
We make the 3+1 identifications 



u° = 7, 
u = 7v, 



(3.3) 



and we note that the spin vector in its (non-inertial) rest 
frame is 



= A»s» 



(3.4) 



where A M v is the Lorentz transformation to the rest 
frame of the particle 



a" , = <>" , + ( 7 1 (7 -7) V 

-7V w 2 ' vv 



(3.5) 



V = U> X V. 



(3.10) 



In the case of uniform circular motion, Eq. (|3.9II becomes 

^=-(7-l)wxso = n r xs , (3.11) 
dt 

where SIt is the Thomas frequency. 



IV. STRING WITH ONE FIXED AND ONE 
MASSIVE END 

A. Spinless quark 

A string with one fixed end and a massive quark on 
the other end is described by an action that is the sum 
of the free massive point particle action and a free string 
action, which we take in Polyakov [22( form, 



5' 



T 



dr I daV-hh^daX^dtX^ 
2 ./ Jo 



— m 



dr v— x 2 . 



(4.1) 



Here X fl (a, r) are the coordinates of the string world- 
sheet parametrized by r and a, h a ^ is the metric on the 
string worldsheet with h = det(/i a &), x^(t) are the co- 
ordinates of the quark worldline, T is the string tension, 
and m is the quark mass. We use small latin letters for 
worldsheet tensor indices. 

We require that the string end at a = is fixed at the 
origin, X(0, r) = 0. To make this an interacting theory, 
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we must impose the condition that the end at a = 1 ends 
on the quark: 



X^(1,t) = x*(t). 



(4.2) 



The variation of the action under variations that pre- 
serve the endpoint conditions, 



6X^(0, t) = 0, 
6X^(1, t) = fe^(r), 



(4.3) 
(4.4) 



SS = dr 



V—i 2 ■ 



T J da dr V^hh ab d a (SX^)d b X„, 

d_ / rnifj, y 
cr=i dr V Y 7 — i 2 / - 

J dadT6X»d a p«, 



drSx^ 



(4.5) 



after an integration by parts. Here we have used the no- 
tation p^ for the current density of spacetime momentum 
on the worldsheet, 



5S/6(d a X») 



T^ 



-hh ab d b X^. 



(4.6) 



We see in Eq. I|4.5|l that the force that moves the quark 
arises from the boundary condition 14.2(1 . The key idea of 
our work is to make a parallel construction with fcrmionic 
variables in the case of a spinning quark. In our con- 
struction, the motion of the quark's spin comes about 
as a result of introducing new fermionic variables on the 
string and the boundary conditions imposed upon them. 



B. Spinning quark 

In this section we make an ansatz for the form of the 
action. In order to have pure Thomas precession, we need 
an action for the fermionic variables £ M whose variation 
has the form 



5 Sot / drS^ ii" +iu"(u-0 



(4.7) 



so that we obtain Eq. I|2.39[l . the condition necessary for 
Thomas precession. 

The term iSt; ■ ii(u ■ £) looks like — i8£ ■ jra~ —iS£, ■ 
p 1 ^/™, if wc use the equations of motion u ■ £ = —£5 and 
make the identification mii^ = . 

We can obtain such a boundary variation by introduc- 
ing worldsheet fermionic variables E fJ, (a, r) and E^er, r) 
whose boundary conditions are 

S"(l,r) = e(r), 

£ 5 (1, t) = &(r), (4.8) 
and then replacing d a X^ in the string action (|4.1|l by 

W a = d a X»-a-d a E»E 5 -P-E»d a E 5 . (4.9) 
m m 



We will fix the parameters a and /3 by requiring con- 
sistency of the equations of motion and pure Thomas 
precession of the spin. 

In analogy to the spinless case, we take our action to 
be the sum of the free Berezin-Marinov [l2| action (|2.16() 
for the particle and a Polyakov action modified by the 
replacement of d a X^ by 11^ defined in Eq. I|4.9|) : 



5' 



dr 



T 



— dr da 



-hh ab WIl 



b /Y 



(4.10) 



C. Equations of motion 

Under variations of the X^ and x^ L that obey the 
boundary conditions Eqs. (|4.3|) . we find the variation of 
the action to be 



SS = dr 



— + 1\ 



\/—x'- 



T J dadTV^hh ab d a (6X")U btl , 

_ d_ / rax^ j\ P ^ V \ 
<7=l dr \y/-x 2 \J-x 2 > 



= \ drSx^ 



da dr 6 Xi 1 d a p°, 



(4.11) 



where we have again used the notation p° for the current 
density of spacetime momentum on the worldsheet, which 
in this case is 

p^ = SS/S(d a X^)=-TV Z hh ab Il blx . (4.12) 

The vanishing of the variation SS leads to equations of 
motion for the quark and the string, 



^ = F, = -TV^hh lb n b , 
= d a (V^hh ab n b ^, 

where the quark's momentum is given by 

P f v 



= mu^ + iX 



\/—x 2 



(4.13) 
(4.14) 

(4.15) 



with the usual projector, — r\^ v + u^u v . 

Under variations of the fermionic variables E^(a, r) 
£5(0-, t), ? m ( t ) and &>(t), obeying 



5S"(l,r) = 5e(r), 
SE 5 (1,t) = S£ 5 (r), 



(4.16) 



that preserve the boundary conditions l|4.8|l . we find the 
variation of the action to be 



6 



SS = dr 



iX S^ 1 Up + iX 



m 



da dr 



a(d a 6E" E 5 + d a E» SE 5 ) + p(6Z» d a E 5 + Z» d a SE 5 



dr 



[iX Ufl - if M + a-J-£ 5 F M ) 6^ + (iX - i& - F„ ) 



(It 



(a£S"S 5 +/3H^H 5 ) 



cr = 



— / da dr 

to 



6E 5 



(4.17) 



Using the notation of Eq. (|4.13[> . and the equation of 
motion l|4.14|l . we find the equations of motion 



C = Xu>* + a& — , (4.18) 
to 

6 = A-/3£ M — (4.19) 
m 

= (a-/3)\/=ft/i o6 n o „abE 5 , (4.20) 

= (a- l3)V^hh ab H a ^d b E^ (4.21) 



These last two equations of motion, Eqs. I|4.20|l and 
(I4.21|) . would be automatically satisfied if a = /3. 

The equation of motion for the metric h ab yields the 
vanishing of the stress-energy tensor, also known as the 
Virasoro constraint, 

T ab = n^U b ^-h ab h cd U^U d ^0. (4.22) 

Variation of the multiplier A yields the equation of mo- 
tion, 

(4.23) 

that becomes the Dirac equation constraint in canonical 
language 

Ppt" + mfc « 0. (4.24) 
The Klein-Gordon mass-shell condition, 

l(p 2 + m 2 )«0, (4.25) 

arises directly from squaring the momentum l|4.15|l . 
Equation (|4.25() can also be found by taking the Dirac 
bracket of the constraint 14.2411 with itself, as in 
Eq. pTEty . 

We also need boundary conditions on the string 
fermionic variables at the fixed end, S M (0, r) and Ss(0, r), 
in order to make the second integral in Eq. 14.17( 1 vanish. 
We cannot impose = TV— hh u H b ^\ a= o because that 
is the force on the fixed end, which cannot vanish. The 
correct boundary conditions are Dirichlet, of which the 
simplest are 

S"(0,r) = 0, (4.26) 
S 5 (0,t) = 0. (4.27) 



V. DETERMINATION OF a AND /3 

A. Conservation of the Dirac equation constraint 

We begin by looking at the equations of motion in A = 
gauge in order to make the ideas clearer. With A = 0, 
Eq. 1 (4.13( 1 becomes 

u" = — . (5.1) 

TO 

Using this, we simplify Eq. (|4.18jl and Eq. 1(4.19(1 to 

e = a&u", (5.2) 
is = (5.3) 

The equation of motion 1(4.23(1 that leads to the Dirac 
equation, 

u%, +& = (>, (5.4) 
must be constant in time for consistency. We find 

= u ■ £ + a^u ■ u — /?£ • ii 

= (l-/3)u-£ = 0. (5.5) 

Thus, for consistency we must have (3 = 1- 

In a general gauge with A ^ 0, we obtain a similar 
result: 

d n ■ F F ■ f 

JL(p^ + mt 5 ) = a £ 5 ^- + (l-/3)^ 
dr mm 

= (1-/3)^-1=0, (5.6) 

TO 

as long as p ■ F — 0, which is required for the consistency 
of the mass-shell relation (|4.25() . We take up this issue 
at the end of this section. 

B. Thomas Precession 

Using the equation of motion (|4.23|) in 15.211 , we find 
= -mt"(if£) +Xu". (5.7) 
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The analysis of Sec. HT1 showed that it was necessary for 
Eq. 12.39(1 to hold in order have pure Thomas precession. 
Comparing Eq. (|5.7(l to ((2.39(1 . we find it necessary that 
a = 1 in order to have pure Thomas precession. 



C. Consistent action and boundary conditions 

Because a = (3 = 1 from the consistency and pure 
Thomas precession requirements, the string variable 11^ 
is a total derivative, 



K = daX^ 



with 



X* = x>- 



(5.8) 



(5.9) 



Remarkably, this combination is also the key to simpli- 
fying potential interactions of two fermions |23| . 

The consistent action for a QCD string with a spinning 
quark on one end that undergoes pure Thomas precession 
can be written using 1)5.9(1 as 



5' 



dr 



J o d(J J dTy/^hh ab d a X' 1 d b X tl . 



(5.10) 



Because we have a = f3, the equations of motion ((4.20(1 
and 1(4.21(1 are automatically satisfied and the boundary 
conditions on S M and S5 at the fixed end can be relaxed 
slightly; 



S"(0,t) Sb(0,t)=0. 



(5.11) 



D. Conservation of the mass-shell constraint 

We used the condition p ■ F = in Eq. 1(5.6(1 . This 
condition is also necessary for the conservation of the 
mass-shell relation ((4.25I) ; 



1 d , o 
= - — (p 2 + m 2 ) = p ■ F. 



(5.12) 



We show that Eq. ((5.12(1 follows from the equations of 
motion of the full action 1(5.10(1 . To begin, we use the 
equations of motion 1(4.13(1 . 1(4.18(1 . and ((4.19() and the 
expression (I4.15|) for the quark's momentum to calculate 
the boundary value 



k 

m 



i(\u^ + —F^)c 5 
\ m J 



if A- 



m 
i 



V^i 2 p" + — ■ F. 



(5.13) 



Using the nilpotency of £ • F and the Virasoro constraint 
KT^i . we find 



p-F 



F, 



my—hh 



ib 



[m\/—hh 



(7=1 



lb 



= 

= 0. 



-h ob h cd d c X»d d X, 



(5.14) 



We have used h lb hoi, = Sq = in the last line. If we do 
the same analysis keeping a and [3 arbitrary, after a bit 
of algebra we find 



p-F = i(l-0)-jL=(£-F), 
again showing the necessity of having /3 = 1 . 



(5.15) 



VI. FERMIONIC GAUGE INVARIANCE 

The string portion of the action 1(5.10(1 has two 
fermionic constraints, 



= n M - -p^-s « 0, 

m 



$ 5 = n 5 + -i>5" « 0, 



(6.1) 
(6.2) 



where ILj, II5, and are the momenta conjugate to 
S^ 1 , S5, and respectively. The fermionic constraints 
together with the Virasoro constraints ((4.22|) are all first- 
class. It is easy to compute the Poisson brackets 

$ 5 } = {$ 5 , $ 5 } = 0. (6.3) 

Because the stress tensor ((4.22(1 is traceless, there are only 
two independent Virasoro constraints, which we may take 
in the form 

L ± = X -(P±TX'f . (6.4) 

After a bit of algebra, we find the Poisson brackets 

{L ± (a),L ± (g)} = T (L±(a) + L±(g)) S'(a — g), 

{L ± (a),L T (g)} = 0, 

{L ± (a),%(g)} = 0, 

{L±(a),$ 5 (g)} = 0. (6.5) 

By acting in combination on the fields A through Poisson 
brackets, 



6 H A = {A, HV^+H&s}, 



(6.6) 
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the constraints l|6.1(l generate the following fermionic 
gauge invariance of the action 



6 H X" = -(H»~ 5 + ~»H 5 

777 V 



(6.7) 



Here if = (a, r) and ffg = H5 (a, r) are Grassmann- 
valued functions on the string worldsheet. This is not an 
invariance of the particle action, so the gauge parameters 
H must vanish at the boundary. Obviously X 11 is gauge 
invariant, 



ShX* = 0, 



(6.8) 



so the string action i|5.10[) is invariant as well. Because 
we have as many first-class constraints as fermionic vari- 
ables, there are no dynamical fermionic degrees of free- 
dom on the string; except for their values at the bound- 
ary, they are pure gauge. 

Just as for a free Dirac particle action, the particle 
piece of the action (|4.1U|I has a local supersymmetry gen- 
erated by the Dirac constraint 

(j) = p ■ £ + 777^5 ~ 0. (6.9) 

The gauge variation of x^ is 

S n x^ = {x»,ir)4>} D = irj^. (6.10) 
The gauge variations of the other variables are 



5 V P» = 0, 
^£5 = -rym. 



(6.11) 



The Lagrange multiplier fields A and e have gauge varia- 
tions 



5 n X 



-77777, 
27.A77 



(6.12) 



VII. ENERGY AND ANGULAR MOMENTUM 

The numerical quantization of the relativistic flux tube 
model starts from the conserved quantities of the system. 
The canonical form of these quantities for the string with 
a spinning quark are only slightly different from those for 
the spinless case. In this section we calculate the four- 
momentum and the angular momentum of our system. 

The action, Eq. JOB, is invariant under infinitesimal 
translations and Lorentz transformations, 



<fcc"(T) 

<DC' 1 (ct,t) 

se(r) 



^vr(r), 



Noether's theorem guarantees the existence of conserved 
total momentum and conserved total angular momen- 
tum J^v, which can be computed from the vanishing 
change in the action, 

5S = aPKPt + ^AJ^ = 0, (7.2) 

assuming the use of the equations of motion. From 
Eq. H7.2fl . we find explicitly 

■P" = p"+ f daP" (7.3) 
Jo 

Jp,v = x\)jp v ] ~ C[fj,^u] + J da (XfrP,,] - S[ M n 

= x^jp v ] - i^iu + I do X {li P u] , (7.4) 
Jo 

where X v = X v - ^S V S 5 . In the last line of Eq. fTljl 
we have used the constraints 12. 2[) and JOJ- 



VIII. DISCUSSION 

Starting from the requirement that an action for a spin- 
ning quark on the end of a string should lead to pure 
Thomas precession of the quark spin, we have shown 
how to construct a consistent action for a massive spin- 
ning quark on the end of a QCD string whose other end 
is fixed. To do so, we introduced additional fermionic 
variables on the QCD string itself and required that the 
equations of motion of the quark's fermionic variables 
arise from boundary conditions on the string. The two 
parameters we introduced, a and /?, are fixed to unity by 
the requirements of pure Thomas precession and consis- 
tency of the equations of motion, respectively. 

Our action has a consistent fermionic symmetry that 
allows us to gauge away the stringy fermionic degrees of 
freedom. Other authors 0> HE HH have argued for the 
existence of a supersymmetry on the QCD string world- 
sheet when there are spinning quarks on the boundary. 
It would be interesting to know if our fermionic symme- 
try is that supersymmetry. As one piece of evidence, we 
can make contact with the Wilson super-loop approach 
|25ll26U27| by making a change of variables in our action, 
Eq. H5.10fl . When we make the variable change 



rn 
i 



^ '5 ; 



(8.1) 



in the string and quark actions respectively, we are led 
to the Polyakov bosonic string action plus a quark action 



S q = J dr (^—m\J~^i 



x 2 - ivPii 1 ' £u.£,v + 



(8.2) 



(7.1) 



where we have neglected to write the kinetic terms for 
the fermions, the Dirac constraint, some higher-order 
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fermionic pieces, and a total time derivative. The sec- 
ond term yields the interaction of the quark's spin with 
the string worldsheet. Unfortunately, this action is not 
simple and does not seem to have a reasonable canon- 
ical formulation, unlike Eq. (|5.10() . We also note that 
the second term of Eq. 1|8.2[) is similar to one added by 
the Martemyanov and Shchepkin |15| . though there are 
additional terms in (|8.2[1 not present in their action. 

We have not considered the case of a quark at each end, 
however the generalization is immediate. We introduce 
a set of Grassmann variables for each quark, £ 5i , 
with i = 1,2. In this case, however, the variables for 
one quark commute with those of the other [2j|, just as 
the gamma matrices of two different fermions commute. 
We also introduce a set of worldsheet fermionic variables 
for each quark and make the generalization in the string 
action 



(8.3) 



i=l,2 



Though the S variables of each quark are Grassmann- 
valued, the 5 variables of one quark should commute 
with the S variables of the other. 

We have not considered the case of a massless quark, 
which appears to be somewhat problematic in our for- 
malism. On the other hand, there is no problem treating 
very light, but still massive, quarks in this formalism. 

Although we have partially analyzed our action in the 
general case of a 7^ /3, we have not pursued the analysis 
with a 1 because we are most interested in the phe- 
nomenologically relevant a = /? = 1 case. In the more 
general case, the fermionic constraints on the worldsheet 
are not all first-class and some of the fermionic variables 
on the worldsheet may become dynamical, though addi- 
tional terms in the action may be necessary to preserve 
the first-class nature of the Virasoro constraints. 

We hope to present soon a numerical quantization of 
this action. 



The variation of H under a change of a fermionic variable 
such as S^ a is 



d R H 



(A2) 



where d R H/d^ a denotes the derivative from the right. 
We could equally well have used 



SH = 6? 



,d L H 



(A3) 



since it has the same value. Variation of the action leads 
to the usual canonical equations of motion, 
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Pi 



dH 

dpi 



dH 








d R H 


d L H 


dir a 


dir a 


d R H 


d L H 







(A4) 



The last two relations of Eq. (|A4|) follow because H is an 
even Grassmann parity function. These relations can be 
succinctly summarized by the introduction of a Poisson 
bracket 



= {z,H}, 



(A5) 



where z is any of q l , pi , £ Q , or n a and the Poisson bracket 
of any two functions A and B is 



{AB} = ]T 



dA dB dA dB 

dq l dpi dpi dq l 



fd R Ad L B d R A d L B\ , A „ S 

(A6) 
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APPENDIX A: PSEUDOCLASSICAL 
MECHANICAL CONVENTIONS 

We take a canonical form for an action to have the 
velocities to the right of the momenta, 



S = Ldr= dr 



p l q l + TT a i a - H{q,j>,ir,£) 

"(Al) 

where q l and pi are bosonic variables and £ Q and 7r Q are 
fermionic variables and H is a Grassmann even function. 



When a system has second-class constraints and one 
wishes to set them strongly to zero, consistency requires 
that the Poisson brackets of the system be modified so 
that the Poisson bracket of any second-class constraint 
with any other phase space function is identically zero. 
This modified Poisson bracket is called a Dirac bracket. 

The simplest example is illuminating, though artificial. 
We imagine a dynamical system in which there are 2N 
phase space variables and two second-class constraints, 
qN w and pat « 0. These constraints hold for all time so 
qN and pn are irrelevant variables; no physical quantity 
should depend upon them. The correct procedure is to 
ignore these variables and replace the Poisson bracket, 



N 



dA dB 

dq l dpi 



dA dB 

dpi dq 1 



(Bl) 
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by the Dirac bracket 



N-l 



»=i 



dA_dB_ _ dA dB 

dq i dpi dpi dq i 



(B2) 



The Dirac bracket (|B2ll of q^ or p^ with any other phase 
space function will obviously vanish. 

For a more complicated system with second-class con- 
straints Xi ~ 0j the Dirac bracket is less obvious. The 
matrix of Poisson brackets of the second-class constraints 

{ X i,Xj} = A ij , (B3) 

has non- vanishing determinant, and is therefore invert- 



ible. We denote the matrix inverse to Ajj by A 1 -? , and 
define the Dirac bracket of any two functions A and B as 

{A, B} D = {A, B} - {A, x% }A^{ Xj ,B}. (B4) 

The desired property now follows, 

{A, X k}D = {A, X k}-{A,Xi}^ ij {Xj,Xk} 
= {A, Xk } -{A, X i} A«A jk 
= {A, X k}-{A, Xt }S t k = Q. (B5) 

We note that some authors use {A, B}* to denote the 
Dirac bracket. 



[1] M. B. Green, J. H. Schwarz, and E. Witten, "Superstring 

Theory. Vol. 1: Introduction," (Cambridge Univ. Press, 

Cambridge, 1987). 
[2] M. Ida, Prog. Theor. Phys. 59, 1661 (1978). 
[3] D. LaCourse and M. G. Olsson, Phys. Rev. D 39, 2751 

(1989); M. G. Olsson and S. Veseli, Phys. Rev. D 51, 

3578 (1995). 

[4] W. Buchmiiller, Phys. Lett. B 112, 479 (1982). 

[5] A. B. Henriques, B. H. Kellett and R. G. Moorhouse, 

Phys. Lett. B 64, 85 (1976); H. J. Schnitzer, ibid. 65, 

239 (1976); 69, 477 (1977); Phys. Rev. D 18, 3482 (1978); 

L. H. Chan, Phys. Lett. B 71, 422 (1977). 
[6] N. Isgur, Phys. Rev. D 57, 4041 (1998). 
[7] E. Eichten and F. Feinberg, Phys. Rev. D 23, 2724 

(1981); D. Gromes, Z. Phys. C 26, 401 (1984). 
[8] G. S. Bali, K. Schilling, and A. Wachter, Phys. Rev. D 

56, 2566 (1997). 
[9] J. L. Martin, Proc. Roy. Soc. A 251, 536 (1959); 

R. Casalbuoni, Nuovo Cim. A 33, 115 (1976); ibid., 389 

(1976) . 

[10] P. A. M. Dirac, Lectures on Quantum Mechanics, 
Belfer Graduate School of Science (Yeshiva Univ., New 
York, 1964); K. Sundermeyer, Constrained Dynamics, 
(Springer, Berlin, 1982). 

[11] M. Henneaux and C. Teitelboim, Quantization of Gauge 
Systems, (Princeton Univ. Press, Princeton, 1992). 

[12] F. A. Berezin and M. S. Marinov, Annals Phys. 104, 336 

(1977) ; JETP Lett. 21, 678 (1975). 



[13] R. Casalbuoni, Phys. Lett. B 62, 49 (1976); A. Barducci, 
R. Casalbuoni, and L. Lusanna, Nuovo Cim. A 35, 377 
(1976). 

[14] L. Brink, S. Deser, B. Zumino, P. Di Vecchia, and 
P. S. Howe, Phys. Lett. B 64, 435 (1976); L. Brink, P. Di 
Vecchia, and P. S. Howe, Nucl. Phys. B 118, 76 (1977). 

[15] B. V. Martemyanov and M. G. Shchepkin, Z. Phys. C 
48, 647 (1990). 

[16] J. Frenkel, Z. Phys. 37, 243 (1926). 

[17] L. H. Thomas, Phil. Mag. 3, 1 (1927). 

[18] V. Bargmann, L. Michel, and V. L. Telegdi, Phys. Rev. 
Lett. 2, 435 (1959). 

[19] F. Ravndal, Phys. Rev. D 21, 2823 (1980). 

[20] C. A. Galvao and C. Teitelboim, J. Math. Phys. 21, 1863 
(1980). 

[21] J. Gomis, K. Rafanelli and M. Novell, Phys. Rev. D 34, 
2298 (1986). 

[22] A. M. Polyakov, Phys. Lett. B 103, 207 (1981). 

[23] P. Van Alstine and H. Crater, J. Math. Phys. 23, 1697 

(1982); Annals Phys. 148, 57 (1983); Phys. Rev. Lett. 

53, 1527 (1984); Phys. Rev. D 33, 1037 (1986). 
[24] P. Horava, Nucl. Phys. B 463, 238 (1996); JHEP 9901, 

016 (1999). 

[25] A. M. P olyakov, Nucl. Ph ys. Proc. Suppl. 68, 1 (1998). 
[26] V. Vyas, hep-th/0010166 

[27] A. Migdal, Prog. Theor. Phys. Suppl. 131, 269 (1998). 



